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Abstract 

The interlayer pair-tunneling (ILPT) mechanism for high-Tc superconductivity is 
able to predict the dependence of the (optimal) critical temperature on the num- 
ber of layers n within an homologous series of layered cuprate oxides. We generalize 
the mean-field procedure employed to evaluate within an extended in-plane Hub- 
bard model in presence of ILPT, developed for a bilayer complex (n = 2), to the 
case of n = 3, 4 inequivalent superconducting layers. As a function of doping, we 
show how a nonuniform hole-content distribution among different layers affects Tc- 
In particular, depending on doping, the onset of superconductivity may be ruled by 
inner or outer layers. The latter result may be related to recent experimental data 
of Tc as a function of pressure in Tl- and Bi-based layered superconductors. 
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1 Introduction 



Among the features that have been recognized to characterize most high-Tc 
cuprate oxides are the dependences of the critical temperature Tc on hole- 
doping 6, and on the number n of Cu02 layers in Bi-, T1-, and Hg-based 
multilayered compounds. While Tc = Tc{6) is usually well represented by 
a bell-shaped curve, Tc is seen to increase monotonically with increasing n, 
tending to a different saturation value for each homologous series for n ^ 1. 
While the former fact strongly supports the idea of a bidimensional pairing 
mechanism, the latter has suggested [1] that coherent tunneling of supercon- 
ducting pairs between adjacent Cu02 layers may considerably enhance Tc, and 
provide the in-plane order parameter Ak with an unconventional anisotropy 
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in the Brillouin zone, its symmetry being determined solely by the nature and 
symmetry of the in-plane pairing [2]. A basic requirement is that the uncon- 
ventional properties of the normal state suppress or greatly reduce coherent 
single-particle tunneling. This may be ensured by Anderson's orthogonality 
catastrophe or by spin-charge separation [1,2]. 

In order to single out the relevance of interlayer against purely bidimensional 
superconducting mechanisms, much theoretical as well as experimental effort 
has been devoted to the study of the two hmiting cases of infinite-layered 
(n = oo) compounds, and of single and bilayered compounds, with almost 
isolated layers. Another means to investigate the competition between 2D 
and interlayer effects is naturally provided by multilaycrcd cuprates, with an 
intermediate number (n = 2, 3) of, but inequivalent layers. 

Within the ILPT mechanism, one source of inequivalence among layers would 
be provided by the mechanism itself. In fact, a given Cu02 layer may be 
coupled through ILPT to either one or two adjacent layers, depending on its 
position within a single unit cell. Therefore, one has to distinguish between 
inner and outer layers. This effect alone is able to account for the observed 
dependence of Tc on n [1] . Moreover, two layers are made inequivalent by their 
different position with respect to the 'reservoir' blocks at the edges of the unit 
cell, which induces a nonuniform distribution of hole-content among them. As 
said above, this is expected to affect bulk superconducting properties such as 
T 

c- 

In the following, we shall address the issue of the competition between these 
two sources of inequivalence in the case of a multilayered complex with a 
slightly nonuniform hole-content distribution. 



2 The model 

We consider the following model Hamiltonian for tightly bound interacting 
fermions in an n-layered complex [3]: 

— Sk'-ko-'-ko- Z_/ '^kk''-kT'--kJ,'--k't'^k'|- \^) 
kat Vk'W 

Here, o^^ (c^^) creates (destroys) a fermion on layer ^ {£ — l,...n), with 
spin projection a I along a specified direction, wave- vector k belonging to 
the first Brillouin zone (IBZ) of a two-dimensional (2D) square lattice, and 
band dispersion = — f^^, measured relative to the chemical potential 
H^. At variance with Ref. [3], here we explicitly assume that may take 
different values depending on the layer index i, thus accounting for different 
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hole content (band filling) nf^ on incquivalent layers [4] . Furthermore, on each 
layer we assume the 2D tight-binding dispersion relation (a being the lattice 
step): 



where at least nearest- neighbours (t = 0.25 eV) as well as next-nearest neigh- 
bours it' /t = 0.45) hopping has to be retained, in order to reproduce the most 
relevant properties common to the majority of the cuprate compounds [5] . 

In the Hamiltonian Eq. (1), interaction is restricted in the singlet channel only 
through: 



where 9{r) is the usual Heaviside step-function. The interaction term is thus 
made of an in-plane contribution t/kk', which provides for Cooper pairing 
within each layer, and an effective contribution, arising from coherent pair- 
tunneling, here restricted between adjacent layers only. The ILPT matrix el- 
ement Tj(k) describes a second-order perturbation in the hopping matrix 
element t_L(k) orthogonal to the CuO layers. Therefore, T/(k) = t\{k)/t, 
where [2] t±{k) depends on k locally in the IBZ as t±{k) = t±[cos{kxa) — 
cos(/Cj^a)]^/4, as recently confirmed by detailed band-structure calculations [6], 
with t± = 0.1—0.15 eV. 

We shall not attempt at specifying the microscopic origin of the in-plane pair- 
ing mechanism. However, the potential f/kk' has to be invariant under the 
symmetry transformations of the underlying crystal point group. Therefore, 
it may be expanded as a bilinear combination of basis functions for the ir- 
reducible representations of such group, viz. C^y for the 2D square lattice. 
Such a series truncates after a finite number of terms, in the case of a finite- 
ranged potential. In Ref. [7] we have considered the competition among differ- 
ent symmetries arising from such a representation, in the presence of ILPT, as 
a function of band filling. Since we shall be mainly interested in the optimal 
filling region, where d-wave has been proved to prevail [7] , in agreement with 
most experimental results, we may take ?7kk' = '^3fl'3(k)fi'3(k')- Here, following 
the notation of Ref. [7], 5'3(k) = ^[cos{kxa) — cos(A;j^a)], and A3 < is a phe- 
nomenological constant, related to in- plane inter-site coupling [7]. We remark 
that, in any case, no symmetry mixing would occur at exactly the critical 
temperature [7]. 

Standard mean-field (MF) treatment of Eq. (1) yields the BCS-like gap equa- 
tion: 



£k = —2t[cos{kxa) + cos(A;j,a)] -|- At' cos{kxa) cos{kya), 



(2) 



V^' = 4^kk'(^«' - Tj(k)(5kk'(l - - K - ^1), 



(3) 
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^k— '^kk'Ak'^k' 

k'£' 

-I 

k' 



—It. Ur^^'xlK' + TAk) \xi^W^^' + xi-'K-'] , (4) 



where Xk denotes the pair susceptibihty on layer i. In the case of a multilay- 
ered complex, one usually introduces two more auxiliary gap functions A^, 
identically vanishing, on the fictitious layers i = and £ = n + 1 [3]. Eq. (4), 
supplemented by these conditions, explicitly displays two sources of inequiva- 
lence between layers. One source is the different value that /i^ may assume on 
different layers. This is employed to describe a nonuniform hole-content dis- 
tribution among layers. On the other hand, superconductivity in inner layers 
(1 < £ < n) is enhanced by ILPT with two adjacent layers, whereas pairs in 
layers i — 1 and i — n, i.e. at the bottom and at the top of the n-layered 
stack, respectively, can tunnel coherently only towards one adjacent layer. It 
is known that both doping and the ILPT mechanism separately contribute 
to determine the critical temperature [7]. In the following, we shall study the 
interplay of these effects in the case of an n-layered complex. 

In order to account for 'edge effects' [3], the Ansatz A^ = Aksin(£7) is in- 
troduced, thus allowing to decouple the k-space from the ^-dependence of the 
order parameter [8]. The condition A^ = A|^^^ = then yields 7 = Lin- 
earization of Eq. (4) at T = Tc yields the following n nonlinear equations for 
the critical temperature Tc'. 



Here Xk — (2Ck)~^ ^^^'^(|/^c^k) denotes the linearized pair susceptibility for 
layer £. The renormalization factor, due to the ILPT mechanism [2], here 
contains the 'averaged' pair susceptibility 



Xk 



sm 



'Jn_ 
n + 1 



n -1 r 



, (6) 



which reduces to the expression obtained in Ref. [3], when /i is the same for 
all layers. 

For each i — 1, . . .n, Eq. (5) imphcitly defines the critical temperature 
corresponding to the superconducting instability on layer i. Eqs. (5) are de- 
coupled with respect to £, by virtue of the assumptions above, and because of 
linearization at T^. However, the whole n-layered stack undergoes a metal-to- 
superconducting transition as soon as one of the n layers does, at least within 
this simplified MP approximation. Therefore, the largest of the solutions 
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of each of Eqs. (5) is to be regarded as the true critical temperature Tc- 
One can straightforwardly recognize that the equation 

min[l - Tj{k)xi] = 0, = Tf , (7) 

k 



implicitly defines a lower bound to T^f. In analogy with the bilayer case [7], 
Eq. (7) can be solved analytically for as a function of //^, showing that 
the ILPT mechanism actually enhances superconductivity at all band fillings. 
In the case of uniform hole-content distribution, /i^ = fi, all i, near optimal 
doping Eq. (7) reduces to the known expression T^/Tj = |cos-^^, where 
Tj = t]_/t, in qualitative agreement with the observed dependence of on 
the number of layers n within the majority of layered HTCS [3]. 



3 Numerical results 

Eqs. (5) have been solved numerically for as a function of j/ for n = 3 
and 71 = 4. In these cases one may assume, for reasons of symmetry, that 
jj} = jji^ = = and jj} = jjL^ = n"^ — ji^ = respectively, 

which define the hole-contents n\l^ and n\"'' within inner and outer layers. The 
following values of the parameters have been considered, in order to allow a 
quahtative comparison with the data for Bi- and Tl-based layered cuprates: 
t± = 0.08 eV, A3 = -0.2125 eV [7]. 

Fig. 1 (left) shows as a function of nj^'' and near optimal doping, in the 
case n = 3. In case n = 4, one obtains qualitatively similar results, not shown 
here. According to the above discussion, one can easily identify whether inner 
or outer layers give rise to superconductivity upon varying the hole-content 
distribution among layers. Along the contour n^^^ = nl^\ corresponding to 
imiform hole-content distribution, is the same for all layers, and one recovers 
the typical bell-shaped dependence on doping. On the other hand, a shift from 
the uniform distribution makes either inner or outer layers prevail, depending 
on the actual values of n^*^ and n^jf^ . 

Charge transfer from block 'reservoir' layers is generally agreed to provide the 
superconducting planes with a hole concentration 6, defined as the number of 
holes per unit CUO2. This depends on the crystal structure, as well as on the 
chemical nature and formal valence of its constituents. Reservoirs may also 
act as 'buffers', i.e. they may repel charge excess, therefore confining holes 
within the CUO2 layers. Moreover, oxygen relaxation processes may affect the 
hole-content distribution, due to the relatively high mobility of oxygen defects 
in the cuprate oxides, even at RT [9]. In multilayered cuprate compounds. 
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Fig. 1. Left: Critical temperatures Tc for inner and outer layers in the case of an 
n = 3 layered complex, as functions of the hole-contents. Right: Same quantities, as 
functions of doping S, for nonuniform hole-content distribution, as predicted by the 
point-charge model, with ciht = {m*/me)/e = 0.33. The solid curve corresponds to 
the uniform distribution. 

all these effects are believed to determine a nonuniform hole-content distri- 
bution among different layers, as evidenced experimentally [10]. Hydrostatic 
pressure is known to generally increase the overall hole-content as shown by 
Hall-resistance measurements in a variety of compounds [11], and may influ- 
ence its distribution among layers as well. 

Several models have been proposed in order to estimate the hole-contents n^*^ 
and n^^'' within inner and outer layers, respectively, as a function of doping 
6 [12-14]. Such models generally aim at minimizing the total carrier energy 
in the layered complex, expressed as a sum of band [13,12] or ionization [14] 
energy, plus electrostatic energy. In the point-charge model of Haines and Tal- 
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Ion [13], the charge distribution within a given layer is reahstically described as 
locahzcd on the constituent ions of a lattice unit cell. For several compounds, 
the electrostatic energy can be therefore expressed as a Madelung sum [15]. Al- 
though the point-charge model is relatively system-dependent, its qualitative 
predictions depend mainly on the adimensional parameter ckht = (m*/ me) / e, 
expressing the ratio between the band-effective mass m* in units of the electron 
bare mass rrie, with respect to the background dielectric constant e. 

According to the point-charge model, one finds that in n = 3 layered complexes 
the hole-content distribution is almost uniform, whereas forn = 4 the majority 
of carriers lie in the outer layers, yielding a vanishing especially at low 
doping 6 [13,15]. As a result, more holes are available for the formation of 
pairs, and superconductivity is favoured in the outer layers. On the other hand 
superconducting pairs within inner layers may coherently tunnel towards two 
adjacent layers. Therefore, due to the ILPT mechanism, one expects a larger 
enhancement of within inner than within outer layers. As a function of 
doping 6, the two effects compromise, and a fairly complex scenario arise, as 
shown in Fig. 1 (right). 

Both in the case n = 3 and n = 4 one observes that the onset of supercon- 
ductivity is driven by outer layers at low doping, whereas inner layers rule 
out at higher doping, and particularly around optimal filling, yielding larger 
critical temperatures. Therefore, one may argue that the effect of the ILPT 
mechanism prevails near optimal doping. However, a second crossing may oc- 
casionally occurs at higher dopings, where superconductivity may be again 
dominated by outer layers. 

The behaviour of Tc — Tc{S) displayed in Fig. 1 (right) may eventually account 
for the doping contribution to the dependence of Tc on hydrostatic pressure 
P [16], assuming an approximately linear dependence of 5 on P. The occur- 
rence of 'crossovers' from T^(°^ to T^^*) reminds one of the cusps experimentally 
exhibited by Tc(P) in n = 3 Tl-2223 and n = 4 Tl-2234 [15]. Before a compar- 
ison could be made, however, intrinsic pressure effects should be realistically 
taken into account. In particular, the band structure parameters [16], as well as 
the ILPT amplitude and possibly the in-plane coupling constant, are expected 
to depend intrinsically on pressure. 



4 Summary and concluding remarks 

In conclusion, we have generalized a 2D extended Hubbard model in presence 
of ILPT mechanism for the case of a complex of n inequivalent layers. Different 
layers are made inequivalent by explicitly allowing a nonuniform hole-content 
distribution. The point-charge model has been employed to estimate the hole- 
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contents on inequivalent layers, for the simplest cases n = 3,4. Moreover, the 
ILPT mechanism makes inner layers intrinsically different from outer layers, 
depending on the number (one or two) of adjacent layers. The critical tempera- 
ture has been numerically evaluated at different band-fillings, in the mean-field 
approximation, for the simplified case of a decoupled dependence on k and i 
of the order parameter. The competition between the enhancement due to 
the ILPT mechanism and the effect of a nonuniform hole-content distribution 
has been followed with respect to a variation of the overall hole-content 6. It 
is found that ILPT generally wins out around optimal doping, where inner 
layers are responsible of the superconducting transition. However, a crossover 
may occur at higher dopings, where outer layers may dominate. These results, 
although preliminary, are reminiscent of the nontrivial dependence of Tc on 
hydrostatic pressure observed in Tl-2223 and Tl-2234. 
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